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ABSTRACT 

The UV resonance transitions of neutral argon Ar l A1066 A and of singly ionized phosphorus 
PiiA963 A originated in the damped Lya system (DLA) at Za,hs= 3.3901 towards QSO 
0000-2620 have been detected by means of the UVES spectrograph at the 8.2m ESO KUEYEN 
telescope. So far, this is the first measurement of An and the second of Pii ever performed in 
damped galaxies and in high rcdshift objects. This DLA is well known for having one of the 
lowest metal abundances and dust content, and the lowest fractional abundance of molecular 
hydrogen H2. The mcasmcd Ar abimdancc is [Ar/H] = —1.91 ± 0.09 which is equal to the 
abundances of the other a-chain elements (O, S and Si). The similarity of the Ar abundance 
with the other a-chain elements implies the absence of significant photoionization by either UV 
background or stellar sources along the sightline throughout the damped Lya system. 

Both log(Ar/0) and log(Ar/S) ratios are found close to those measured in the extragalactic 
H II regions and in blue compact galaxies where O is more abundant by at least one order of 
magnitude . This strengthens the universality of the Ar/0 and Ar/S ratios and lends support 
to the existence of a universal IMF. 

The abundance of the non-refractory element phosphorus [P/H] = —2.31 ± 0.10 confirms 
the low amount of chemical evolution in the DLA. This is the measurement of P in the most 
metal-poor material and shows a subsolar [P/Fe] = -0.27 value. The measured ratios [P/Si] = 
-0.40 ± 0.13 and [P/S] = -0.33 ± 0.13 provide evidence for a mild odd-even effect. 

Finally, a stringent upper limit to the population of the ^Pi level in the ground state of O I 
is derived, which provides a lower limit to the physical dimensions of the Zabs = 3.3901 system 
of L > 7 pc. 

Subject headings: cosmology: observations — galaxies: abundances — galaxies: evolution — 
quasars: absorption lines 
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1. Introduction 

Damped Lya (DLA) systems are absorption systems observed in quasar spectra characterized by large 
column densities (A''(Hl)> 2x 10^° cm~^). They span the whole epoch of galaxy formation and are generally 
considered the progenitors of the present day galaxies. DLA systems are of particular interest because 
it is possible to measure a variety of chemical elements in their interstellar gas-phase with unprecedented 
precision (Lu et al. 1996; Prochaska & Wolfe 1999; Pettini et al. 1997, 1999 and 2000). However, the 
observed abundance patterns in the DLAs proved difficult to account for, due to the possibility that 
a significant fraction of some elements contribute to dust grain formation thus depleting the gas phase 
abundances (Kulkarni et al. 1997; Vladilo 1998). At face value the observed Si, S, Cr, Al, Ni abundances 
relative to Fe are similar to those observed in the old stellar population of the Galaxy, with the notable 
exception of Zn which is overabundant compared to iron. However, if the [Zn/Fe] anomaly is interpreted 
as a Fe depletion in the interstellar medium, the observed abundance pattern would be similar to those 
observed in the warm phase of the interstellar medium with solar chemical composition. Zinc has not been 
considered a totally reliable clement by Prochaska et al. (2000), because of its uncertain nucleosynthesis. A 
sample of chemical elements larger than those presently available, and in particular of elements with little 
affinity with dust grains, is required to better address the kind of chemical pattern present in DLAs. 

In this paper we report on the detection of UV lines of the non-depleted species - neutral argon Ar i 
and singly- ionized phosphorus P II - in a damped Lya absorption (DLA) system at ^abs = 3.3901 towards 
the quasar Q 0000-2620. So far, phosphorus has been measured by Outram et al. (1999) in the DLA at 
Za.hs =2.625 towards GB 1759+7539 and an upper limit on phosphorus abundance has been set for the DLA 
at 2;abs= 2.309 towards PHL 957 (Molaro et al. 1998). A possible identification of the An lines in the 
Za.hs = 2.8112 damped Lya system towards PKS 0528-250 was reported by Srianand & Petitjean (1998), 
however, in their spectrum the An lines were considerably contaminated by hydrogen absorption which 
yields a rather uncertain Ar i column density. 

The Ar we derive here is the first one ever obtained in any DLA and one of the few values measured 
by means of unsaturated lines unlike with the interstellar studies in the Galaxy (e.g. Sofia & Jenkins 1998). 
Being a noble gas with very low condensation temperature, Tc = 30 K (Spitzer & Jenkins 1975), neutral 
argon atoms do not form chemical bonds with solid compounds, and the depletion of phosphorus is found 
to be very small in the ISM where variations in refractory elements do occur (e.g. Dufton et al. 1986). 
The dust-to-gas ratio in the DLA at Zabs=3.3901 towards QSO 0000-2620 was found to be less than 0.2% 
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than that in the Milky Way (Levshakov et ah 2000b), and thus the measured [Ar/H] and [P/H] ratios^ at 
^abs — 3.3901 should not be significantly affected by dust depletion, thus providing a unique probe of pure 
nucleosynthetic yields at high redshift. 

The study of phosphorus may eventually allow to clarify the dust formation processes in DLA galaxies. 
Since P and Fe have nearly the same condensation temperature — 1300 K, these species should deplete 
together in equal proportions if dust is produced in dense and hot shells of supernova remnants where 
thermodynamic calculations of condensation are appropriate. A possible disparity between P and Fe 
gas-phase abundances would favour a different model where the dust grains grow in the interstellar medium 
(see Jura & York 1978, and references therein). 

We discuss here the implications of the Ar and P abundances in connection with the chemical evolution 
of such damped galaxies. We underline that for research on the element production history, Ar is a quite 
unique element due to its ability of tracing how strongly the apparent abundances of other species from H i 
regions could be influenced by ionization. 

The paper is organized as follows: in § 2 we briefly describe observations and analyze An and Pii 
proflles. The results on Ar and P abundances are discussed in § 3 (sections 3.1 and 3.2, respectively). In 
Section 3.3, the estimations of the mean gas density and the linear size of the Zabs = 3.3901 system are also 
given. Conclusions arc summarized in § 4. 

2. Observations, data reduction and column densities 

The DLA at Zabs ~ 3.3901 towards QSO 0000-2620 shows high neutral hydrogen column density 
Af(Hi) = -\Q'i-^A\±am cvoT^ (Lu et al. 1996). This system has been studied extensively since it is one of the 
highest redshift DLAs seen in the light of a relatively bright QSO (Levshakov et al. 1992; Molaro et al. 
1996; Prochaska & Wolfe 1999; Molaro et al. 2000, hereafter Paper I; and Levshakov et al. 2000b, hereafter 
Paper II). 

The spectra of QSO 0000-2620 were obtained during the first commissioning of the Ultraviolet- Visual 
Echelle Spectrograph (UVES, see Dekker et al. 2000) at the Nasmyth focus of the ESO 8.2m KUEYEN 
telescope at Paranal, Chile, in October 1999, and have been released for public use. Details of observations 



2[X/H] EE log [iV(X)/iV(H)] - log [iV(X)/iV(H)]o 



- 5 - 



and of data reduction are given in Paper I. We recall here that the resolving power is R = A/AAinstr — 49000, 
which corresponds to a velocity resolution of ~ 6 km s~^ (FWHM). 

Singly ionized phosphorus, which is the dominant state in Hi gas, has three strong UV lines with 
A = 961.041 A, 963.801 A, and 1152.818 A (the corresponding /-values are equal to 0.3489, 1.458, and 
0.2361, according to Morton 1991). We have detected with confidence only the PiiA963 A line which is 
shown in Figure 1 together with the NiiiA1709 A and FeiiA1112 A lines. The PiiA961 A is blended and 
the same occurs for PIIA1152 A as seen from the EMMI spectrum of Savaglio et al. (1997). 

Neutral argon has two strong resonance lines in the ultraviolet region at A = 1048.220 A and 1066.660 
A. Their oscillator strengths are /io48 = 0.257 ±0.013 and /loee = 0.064 ±0.003 respectively from Federman 
et al. (1992). The whole region of the QSO 0000-2620 spectrum containing the two An lines is shown in 

Figure 2 to illustrate both the crowding of the Lya forest and the few windows which were used for the 
continuum placing. The Ar lines detected at Zabs = 3.3901 are zoomed in Figure 3. 

The data used for Ar result from the sum of only 2 spectra instead of the whole 3 since one of them 
turned out to be badly contaminated by a cosmic ray event. The S/N ratio for the An A1066 A line 
therefore is 40. As can be seen from the figure, the Ar i A1066 A line is relatively clean, but A1048 A is 
partly blended with additional absorption. The Ar l line strengths are such that in the Milky Way they tend 
to be saturated even for tenuous lines of sight, but in our DLA system saturation is not an issue because of 
the low metalHcity Z ~ W~'^Zq (Paper I). 

The measurements of metal absorption lines falling in the Lya forest have to cope with the possible 
presence of hydrogen contamination and with a local continuum, drawn by using windows in the Lya forest 
which could be located quite far from the feature under analysis. 

The PiiA963 A line (Figure 1) shows a rather symmetric profile with its blue side reaching the 
continuum, while there could be a small contamination in the red wing. The Ar line looks slightly more 
contaminated with absorption on both sides. The extra absorption around the Ar l line is a continuum-like 
absorption without signature of discrete features, which makes it hard to model with single Hi clouds. 

Our previous analysis of both metal and molecular hydrogen lines in this DLA has shown that the line 
broadening in the ^;abs = 3.3901 system is mainly caused by macroscopic large-scale, rather than thermal, 
motions. All absorption lines from H2 to Znll show symmetric profiles with the same Doppler parameter 
of 6 ~ 10 km s~^ within the uncertainty intervals (Papers I & II). The measured broadening parameter, 
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defined as 6 = \/2 cturb where aturb is tlie one-dimensional gaussian velocity dispersion of ions along the line 
of sight, gives a turbulent velocity dispersion of about 7 km s~^ , a typical value for the Milky Way disk. 
We can therefore assume that the intrinsic profile of the An A1066 A and Pii A963 A lines is the same for 
all neutral atoms and low ionized species from the H i region since they trace the same volume elements. 

To calculate the confidence range for the Doppler parameter wc make use of all metal lines associated 
with the DLA redwards the Lya emission and minimize by varying all the other parameters. The 
calculated Xmin values as a function of b with the la and 2a confidence levels are shown in Figure 4. At 
Xmin = 1-148 the most likely value for b is lO.OOlgJs ^™ confidence ranges are computed following 

the procedure described in Press et al. (1992), here we fix the value of b and minimize by varying all the 
other parameters. 

The Pii A963 A and An A1066 A lines are fitted with a theoretical one-component Voigt profile to 
the observed intensities together with the full set of 6 metal lines, for 5 ionic species, which are located 
redwards the Lya emission, namely the lines Znii A2026.1360 A , Cm A2062.2339 A, A2056.2539 A, Fell 
A1611.2004 A, Sin A1808.0126 A and Niii A1709.600 A, assuming that b is the same for all lines and 
leaving all the elemental column densities as well as the local continuum (AC/C) around the PiiA963 
A and AnA1066 A lines free to vary. We have therefore a total of 7 free parameters and 111 and 105 
points for Pii and An, respectively, thus giving for the former z/ = 104 degrees of freedom, and 98 for 
the latter. All the lines have been previously adjusted in wavelength by using the value of ^abs = 3.390127 
which corresponds to the redshift of the isolated NiiiA1709 A absorption line (Paper I). The calculated 
theoretical profiles were convolved with a gaussian instrumental profile with 6 km s~^ width. 

The best fit for Pii gives N{P ii) = 4.29(±0.39) x lO^^ cm'^, b = 9.96 km s'^ , and AC/C = -0.0778 
for a Xmin = 1-10. The result is shown in Figure 1 by the grey solid curve, whereas points with error bars 
give the normalized intensities, corresponding to a signal-to-noise ratio S/N ~ 15 per pixel. Note that the 
intensities in the bottom panel of Figure 1 are re- normalized in concordance with the local AC/C. The 
dotted curve shows the un-convolved theoretical Pll profile in the range \Av\ < 40 km s~^ . The column 
densities of Zn, Cr, Si, Ni and Fe are consistent within 0.02 dex with those found in Paper I. The errors in 
the Pii column density are of 0.04 dex as can be deduced from Figure 5 (panel a), where the Xmin ^ ^ 
function of the A''(Pii) column density is shown. The variations for AC/C, as well as 6 as a function of 
A'' (Pii) are also shown in panels (b) and (c), respectively. The grey area in panel (b) marks the range of 
possible continuum deviations which correspond to the la confidence interval shown in panel (a). 
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The analysis of Ar relies entirely on ArlA1066 A while ArlA1048 A of which only a portion is 
unblended, is checked for consistency with the line model but it is not used in the fit. The corresponding 
signal-to-noise ratio in this spectral region is approximately equal to 40. Figure 3 shows the best fit for the 
one-component profile (grey continuous curve) and the intrinsic un-convolved profile for only the Ar I A1066 
A line (dotted curve). This figure demonstrates that the best fit obtained for the An A1066 A line provides 
a consistent fit of the red wing of the stronger Ar I A1048 A line. Note that a decrease of 0.1 dex in the Ar i 
column density is sufficient to destroy the perfect match present in the An A1048 A line. A Xmin = 1-138 
is foimd for N{Ari) = 1.044(±0.069) x 10^^ cm^^ b = 10.05 km s^^ and AC/C = -0.0086. Figure 6 shows 
the confidence levels for the An column density, the variations for AC/C, and the Doppler broadening. 

To test the whole procedure we analyzed the FeiiA1112 A line, also shown in Figure 1. This is the 
only Fe II line in the Lya forest which is not strongly blended. We then compared the results with those 
of the FeiiA1611 A, which is placed redwards the QSO Lya emission in a clean spectral region and for 
which a reliable iron column density has already been derived in Paper I. For the Fell All 12 A line we 
used the empirical oscillator strength (/1112 = 0.0062) derived recently by Howk et al. (2000) by means 
of FUSE observations. A Xmin = 1-0299 is found for N{Fen) = 7.095(±0.875) x 10^* cm'^, b = 10.07 
km s~^ and AC/C = —0.066. In this case we used only 6 metal lines from the red part of the quasar 
spectrum excluding FeiiA1611 from the x^-minimization procedure. Thus we had 6 free parameters 
{A^Niii, A"sin, A^Znii, A^Crii, fe, AC/C}, 93 data points, and = 87 degrees of freedom, and calculated Xmin 
as a finction of A^(Fe ii A1112). The calculated confidence intervals are shown in Figure 7. The mean iron 
column density for FeiiA1112 A is 0.02 dex lower than that derived from FeiiAlGll A in Paper I and is 
fully consistent within errors. This shows that we can indeed recover the ion column density accurately even 
though the line is falling in the Lyman a forest with random Lya interlopers and uncertainties associated 
with the continuum drawing. The derived column densities and abundances for Ar and P, together with 
those obtained for other species in previous studies, are reported in Table 1. 

The results presented above differ quite slightly from those obtained from a more commonly used 
analysis in which the continuum is kept fixed while both the broadening and the column density are left 
free to vary. In this approach we obtain for phosphorus iV(Pii) — ioi2-70±o.04 ^.y^-^ ^^^^^ ^ _ 10. 3l]^'^ 
km s^i , while for argon iV(An) = ioi4 02±o.02 ^^-2 ^ ^ 9.95 ± 1 km s^^ when the PiiA963 A and 

AnA1066 A lines were analized separately. The shift of the Pii line is Af = —0.01 km s^^ , and that of 
the An line center is negligible, namely Av = —0.66 km s~^ , which is less than one- half of the pixel width. 



- 8 - 



The fact that both Unes have the same centers and the same 6-parameters which are very similar to 
those of metal lines observed in the low line density portion of the spectrum redwards the Lya emission 
line, shows that the contamination of Ar I and P II lines observed in the Lya forest is not significant in the 
z^t,^ = 3.3901 system towards QSO 0000-2620. 

3. Discussion 

For both argon and phosphorus there are no determinations of the intrinsic abundances in either 
dwarf or giant stars because argon has optical transitions from very high excitation levels, and phosphorus 
shows too faint lines in stellar spectra. P was measured in chemically peculiar stars (Castelli et al. 1997), 
halo horizontal-branch stars by Bonifacio et al. (1995), and sdB by Ohl et al. (2000). However, the 
elemental abundances in these stars are likely modified by atmospheric processes such as diffusion. A 
firm identification of the Pv and Piv has been made by Junkkarinen et al. (1997) in the BAL QSO PG 
0946+301 and a remarkable overabundance of phosphorus with P/C w 6O(P/C)0 was found. 

The lack of such important observational reference complicates the interpretation of the Ar and P 
abundances in the DLA systems. Nevertheless, we may compare the measured [Ar/H] and [P/H] ratios 
with theoretical predictions based on integrating massive star yields over an initial mass function given, 
for instance, by Timmes et al. (1995, hereafter TWW). In this section we consider such issue and other 
physical properties of the Zabs = 3.3901 system. 

3.1. Argon Abundance 

3.1.1. Local stellar ionization 

Argon is an extremely volatile element. It has been detected and studied in the Galactic ISM by means 
of Copernicus, IMAPS and recently by FUSE. The gas-phase (Ar/H) ratio in the Galactic interstellar 
medium is found to be below its solar value with reduction factors varying with the lines of sight. Fitzpatrick 
(1996) argued that abundances of undepleted elements in the ISM such as C, N, and O are below their solar 
values, but this cannot explain the line of sight variations and the considerable amount of the observed Ar 
deficiency, [Ar/H] < — 0.5 dex. Since Ar can be hardly depleted into dust grains, Sofia & Jenkins (1998) 
suggested that the apparent Ar deficiency is caused by a partial ionization of Ar by nearby UV stellar 
radiation. 
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Ar I has an extraordinarily high photoionization cross section for photons with energy higher than the 
IP = 15.76 eV. Its photoionization cross section is about ten times higher than that of Hi over a broad 
range of energies. Therefore, for a cloud that is not thick enough to shield its interior from outside sources 
of ionizing radiation, the UV photons may penetrate the cloud and partially ionize hydrogen, argon and 
other elements. In this case the fraction of Ar i can be lower than the corresponding fraction of H i since the 
absorption cross section of hydrogen diminishes oc for energies hv > 13.6 eV, an effect which may lead 
to an apparent deficiency of argon with respect to hydrogen if one assumes that both elements are mostly 
neutral. 

Reversing the argument presented by Sofia & Jenkins on the An deficiency in the local interstellar 
medium, the coincidence of the Ar abundance with those of O, Si and S (see Table 1) implies the absence 
of significant photoionizing flux from stars along the line of sight throughout the Zabs = 3.3901 damped 
galaxy. Otherwise, the presence of stellar photons with hu > 13.6 eV would decrease the An abundance. 
The same eff'ect is not observable in O i, although neutral oxygen also has a photoionization cross section 
larger than that of H l. Oxygen ionization fractions are coupled with those of hydrogen via resonant charge 
exchange reactions. Therefore, there is observational evidence that in this particular damped system the 
observed ionic species are from their dominant ionization stages in the H I gas and, hence, their abundances 
do not require ionization corrections which may be caused by a hard radiation field in a uniform and 
fluctuating gas density absorbing region, respectively (Viegas 1995, Lcvshakov ct al. 2000a). We would like 
to emphasize once more that our measurement of Ar I rules out the presence of active star formation regions 
and of significant partly ionized gas along the line of sight at ^abs = 3.3901. 

Steidel & Hamilton (1992) identified a galaxy at 2.8" from the QSO 0000-2620 sighthne, which was 
supposed to be responsible for the Zahs = 3.3901 damped system. Giavalisco et al. (1994) found no Lya flux 
in correspondence of this object by means of narrow-band imaging placing a 3a upper limit of 1.2 x 10~^^ 
erg s~^ cm~^. The object identified in Steidel & Hamilton (1992) is probably a high redshift galaxy at 
z 2.8 (Steidel 2000), and therefore the galaxy at z = 3.39 still remains undetected. However, the upper 
limit provided by Giavalisco et al. (1994) in the Lya emission refers to the sky area around the QSO 
and remains valid regardless of the nature of the DLA candidate. The general lack of Lya emission from 
DLAs has been commonly explained by the effect of dust and resonant scattering which can attenuate 
Lya photons (Chariot & Fall 1991, 1993). However, the extremely low dust content in our system with a 
dust-to-gas ratio of less than 10~^ of the Milky Way value, coupled with no detectable effects on the An 
abundance, suggests that most likely the DLA in question is not a site of intense star formation. 



3.1.2. A universal IMF 



Observations of emission-lines in ionized H II regions of spiral and irregular galaxies provide abundances 
for O, S and Ar among the others. None of these elements is expected to be heavily depleted into the solid 
phase as dust (Savage & Sembach 1996). It has already been pointed out that the abundance ratios [S/0] 
and [Ar/0] appear to be universally constant and independent of metallicity over about 2 decades of oxygen 
abundances (e.g. Henry & Worthey 1999). Since O production depends on the progenitor mass while Ar 
and S do not, the implication of this constant ratio is that either the initial mass function is universally 
constant or that the observed elemental ratios are not sensitive to IMF variations. 

Sulfur and argon originate from explosive oxygen burning in Type II supernova events through the chain 
16Q _^ 160 ^ 28gj _^ 4jjg^ ^j^gj^ 28gj 4jjg ^ 32g^ g^^^ 32g _^ 4jjg ^ se^j. Howcvcr, Substantial amounts of 

S and Ar may be manufactured in Type la supernova events (Nomoto et al. 1997). Considering the values 
of Izotov & Thuan (1999) for the extragalactic Hii metal-poor regions, we have log(S/0) = —1.55 ± 0.06 
(rms) and log(Ar/0) = -2.25 ± 0.09 (rms). Similar results, log(S/0) = -1.48 ± 0.06 (rms) and 
log(Ar/0) = —2.32 ± 0.05 (rms), for the same class of objects have been recently obtained by Kniazev et 
al. (2000). 

These observed ratios are comparable with the theoretical predictions of Woosley & Weaver (1995), 
Nomoto et al. (1997), and Samland (1998) in which log (S/0) is « -1.75, -1.9, -1.6, and log(Ar/0) « 
-2.45, -2.7, and 2.4, respectively. The theoretical predictions fall somewhat below the observed range, and 
either an overproduction of O in the models or a significant production of S and Ar in Type la supernovae 
are advocated. 

Our results listed in Table 1 show that the Ar abundance follows closely the other a-chain element 
abundances in the Zahs = 3.3901 system. The absolute oxygen abundance in this DLA is log (0/H) + 12 w 7, 
i.e. one order of magnitude below the value obtained in the most metal-poor galaxies where these elements 
were measured. The measured log (S/0) and log(Ar/0) ratios are of —1.72 ± 0.15 and —2.40 ± 0.15, 
respectively, i.e. very close to the observed ratio in blue compact and spiral galaxies. This strongly favors 
the universality of (S/0) and (Ar/0) ratios. The implication pointed out by Henry & Worthey (1999) that 
either the IMF is universally constant or the stellar mass range which produces these elements is narrow 
enough to make the ratios insensitive to IMF variations is reinforced and further extended to much more 
deficient material detected at ^abs = 3.3901. 
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3.1.3. a-capture versus Fe-peak element ratios 

Most of a particle nuclei are synthesized primarily by Type II supernovae during the early stages of 
the galaxy chemical evolution, while the Fe-group elements are primarily produced by Type la supernovae 
and build up over longer time scales. As a result, in the early phases of chemical enrichment the a-chain 
elements are expected to be relatively more abundant than the iron-peak elements. 

Argon is a typical product of Type II SNe, and like zinc it is not depleted in the interstellar medium. 
This makes the [Ar/Zn] ratio a new and interesting dust-free diagnostic tool for verifying the presence of 
the a-chain clement enhancement. By combining the present abundances with those previously determined 
(see Tabic 1) we have disposal of a set of 4 indicators for the a-chain clement abimdanccs such as O, Si, S, 
Ar, and 4 indicators for the abundances of iron-peak elements such as Cr, Fe, Ni, and Zn. This is a rather 
unique circumstance among the DLA systems studied so far. 

The measured [0,Si,S,Ar/Cr,Fe,Ni,Zn] ratios reported in Table 2 for quantitative comparison provide 
([a-chain/iron-peak]) = 0.18 ±0.03 dex (the weighted mean and dispersion) which becomes 0.12 ±0.02 when 
the slightly deviating Ni is not considered. In the case that Ni resulted a better Fe-peak indicator then 
Fe, Cr and Zn, we would indeed find that [0,Si,S,Ar/Ni] are +0.41, +0.36, +0.29 and +0.36, respectively. 
However, Prochaska & Wolfe (1999) measured Niii A1751.910 A and found [Ni/H] = -2.335 + 0.085, which 
becomes [Ni/H] = —1.952 ± 0.085 when the \oggf values are taken from Fedchak &: Lawler (1999) as we 
did for the Niii A1709.600 A (Paper I). The slightly different abundance obtained from the two Niii lines 
prevents Ni from being considered a safe proxy to the iron-peak abundances. 

An [alpha/iron-peak] ratio of « 0.1 is significantly lower than similar quantity in Galactic stars with 
comparable metallicities, [0,Si,S/Fe] ~ 0.3 — 0.6 dex (McWilliam 1997, and references therein). The 
diff'erence would be even more striking if the high O values of [0/Fe] = 0.7 found in halo stars at [Fe/H] = -2 
by Israelian et al. (1998) and Boesgaard et al. (1999) were confirmed. Implications of low a-chain/iron-peak 
element ratios are discussed in Paper I, and the reader is referred to it for additional details. Here we notice 
that the lack of significant [Ar/Zn] enhancement corroborates our previous results on this system as well as 
the assumption that, on the basis of their relative element abundances, at least some of the DLA galaxies 
seem to undergo a chemical evolution which differs from that of the Milky Way (Centurion et al. 2000). An 
analysis of the zinc abundances in the thick disk stars by Prochaska ct al. (2000) showed that zinc could 
be mildly overabimdant with respect to iron with [Zn/Fe] = 0.09 ± 0.023, a slight difference from the thin 
disk and halo behaviour found by Sneden et al. (1991). However, in our DLA [Zn/Fe] = —0.03 ± 0.06 and 
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[Zn/Cr] = —0.08 ± 0.06, which suggests that in the DLAs zinc tracks Fe and Cr closely. 

3.2. Phosphorus Abundance 

Phosphorus is likely a non-depleted element. Dufton et al. (1986) and Jenkins et al. (1986) found that 
phosphorus is essentially undepleted along sightlines sampling low density neutral gas and is depleted by 
approximately 0.5 dcx in cold clouds. However, P is found strongly depleted in dense star forming molecular 
clouds (Turner et al. 1990). From FUSE observations, Friedman ct al. (2000) obtained [Fc/P] w -0.6 
and -0.7 for the Milky Way and LMC diffuse interstellar clouds, respectively, and argued for a significant 
incorporation of iron into dust grains in both galaxies. The phosphorus abundance provided here is one of 
the few ones available for this element in any astronomical site and the only one in little processed material. 

In our DLA the abundances of the refractory elements Fe and Cr are close to that of the non-refractory 
element Zn, showing the absence of significant dust concentration along the line of sight (Paper I). Such a 
circumstance is extremely rare among DLAs since even in the cases with low element depletion considered 
by Pettini et al. (2000), a difference between Zn and Fe, Cr and Ni has been always observed, with the 
exception of the damped system towards QSO 0454-1-039 where [Zn/Cr] « within the errors. Moreover, 
the extremely low fraction of molecular hydrogen /(H2) = 2 A^(H2)/A^(H)totai — 4 x 10^^ in our DLA 
(Paper II) rules out a significant depletion of P in a variety of molecular compounds. Taking into account 
that the absorbing gas is mostly neutral, as discussed above, we conclude that the observed gas-phase P 
abundance reflects the total amount of P. 

The low abundance of P is close to those of other elements (see Table 1) and supports independently 
the early suggestion that the .Zabs = 3.3901 galaxy has indeed a low metallicity. In the context of the other 
published data this result suggests a possible mild decrease in the damped abundances at high redshift, 
which, however, is not evident when the column density- weighted mean is considered (Prochaska & Wolfe 
2000; Vladilo et al. 2000). 

Theoretically the processes that produce the odd-Z element phosphorus in late stages of stellar 
evolution are not identified clearly. As Margaret Burbidge (1999) pointed out the P synthesis was virtually 
ignored by Burbidge et al. (1957). The only stable isotope is ^^P, which is strongly underproduced in 
the explosive O burning. P is likely produced in the carbon and neon shell burning by neutron capture 
similarly to the parent nuclei ^^Si and ^°Si. No significant P production is found during the explosion 
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phases (Woosley & Weaver 1995). 

Phosphorus abundance, hke that of the other odd-Z elements Na and Al, should be proportional to rj, 
the total neutron excess in a given amount of matter^, which is 77 « 1.5 x 10^'^{(/(q), where « 0.02 
denotes the fraction by mass of matter composed of nuclei with atomic number Z > 6 (Arnett 1971, 1996). 
As the metallicities decrease, the available nuclei which provide the neutron excess decrease too. However, 
in the primordial stellar yields Arnett (1996) found that the odd-Z nuclei are underabundant with reference 
to the even-Z neighbours, but the effect is much smaller than that predicted by simple arguments of linear 
increase of neutron excess with metallicity. Theoretical models by Woosley & Weaver (1995) lead also to 
substantial P production with P yields showing a metallicity dependence which is more pronounced at very 
low metallicities. 

In our DLA the measured ratios [P/Si] = — 0.40 ± 0.13 and [P/S] — — 0.33 ± 0.13 provide evidence 
for a mild odd-even effect, which, however, is less than what expected on the basis of the Woosley & 
Weaver (1995) and Limongi et al. (2000) yields. A similar value for [P/S]=— 0.34± 0.07 has been derived 
by Outram et al. (1999) in the DLA at Zabs =2.625 towards GB 1759+7539, which with [Fe/H]=-1.31 is 
however somewhat more metal rich than the DLA studied here. 

In stellar spectroscopy the element abundances are traditionally discussed as related to iron, [X/Fe], 
and as a function of the iron-to- hydrogen ratio, [Fe/H]. Following this way, in Figure 8 we compare our 
results with theoretical predictions since in our case there is no significant dust depletion in the Zabs = 3.3901 
absorber and, hence, [Fe/H] can be used as a reference ratio. Figure 8 shows the result of such comparison 
between the elements listed in Table 1 (dots with 3 a error bars) and the theoretical prediction (shaded 
area). The width of the shaded area just above a marked element shows the uncertainty range (±0.2 dex) 
for the relative abundance of this element [X/Fe] at fixed metallicity [Fe/H] = - 2.0. The uncertainty of 
theoretical predictions is caused by unknown model parameters as described in TWW. The upper limit 
for the nitrogen-to-iron ratio is set by calculations which includseenhanced convection prescriptions, but 
the lower limit is not defined at low metallicities since the standard TWW set of massive stars does not 
produce significant nitrogen yields. We see from Figure 8 that there exists a general consistency between 
the observed [X/Fe] ratios and the theoretical predictions. The P abundance is slightly suppressed whereas 

the neutron excess is defined as rj = («„ — np)/{nn + '^p), where rin represents the total number of 
neutrons per gram, both free and bound in nuclei, and Tip the corresponding number for protons (Arnett 
1971). 
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Ar behaves like other a-chain elements. However, the same consistency is not found if compared with 
the halo stellar abundance pattern, where a-chain elements are much more enhanced. There are intrinsic 
uncertainties in the M > 30 M© extremely metal-poor massive-star models, which may be important here. 
TWW suggest that their Fe yields may be overestimated by a factor 2, which, if applied in Figure 8, would 
shift all abundances by 0.3 dex along the ordinates. For instance in the Samland (1998) analysis the a-chain 
elements are more enhanced in agreement with what is observed in the Galactic halo stars. 

The P abundance found in our DLA is slightly lower than Fe and the other iron-peak elements showing 
that P was slightly undersolar at low metallicities. Since P is produced almost entirely by massive stars to 
achieve the higher solar value, its production must increase with the metallicity of the gas also to compensate 
the contribution to iron made by Type I SNae. In other words the P yields have to be necessarily metallicity 
dependent. This illustrates well how direct P observations in metal-poor environments can be used to 
constrain the parameter values in the stellar nucleosynthesis theory. 

The obtained results on P and Fc abundances do not seem to agree with the dust formation processes 
in the Galaxy where iron depletion in the diffuse interstellar clouds has been well studied. In Section 1, 
two possible mechanisms of the dust formation processes were outlined, and it was noted that P and Fe 
abundances may be used for choosing one of them. For instance, if the dust grains grow in the interstellar 
medium we could detect a disparity in P and Fe gas-phase abundances as seen in the Galactic diffuse 
clouds. Obviously this is not the case for the Zq\)s — 

3.3901 system. If, however, the dust grains are formed 
in hot shells of supernova remnants - an example of this process has recently been presented by Arendt et 
al. (1999) who revealed Fe protosilicate in the Gas A ejecta - then we should observe a gas-phase depletion 
for both P and Fe, since they have the same condensation temperature. This has not been found in our 
DLA either. Therefore, neither of the two possible mechanisms of dust formation seem to work in the Zahs = 
3.3901 damped Ly-a cloud. Taking into account that both volatile and non-volatile elements do not show 
significant differences in their abundances from the expected quantities, we may conclude that in this DLA 
possible destructive processes are also insufficient to cause some evaporation of the volatile substances. 

3.3. Neutral hydrogen density 

The number density of neutral hydrogen, nui. in the Zabs = 3.3901 system can be estimated by using 
the relative population ratio of the first excited state Ol* ('^Pi) to the ground state Ol ('^P2)- 
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In Hi regions with temperatures ranging from 100 K to 1000 K, the rate coefficients for coUisions 
with hydrogen atoms are larger than those for electron collisions when the electron density, Tie, is small, 
i.e. rie/nni ^ (Bahcall & Wolf 1968). If we also neglect the contributions of proton and molecular 
hydrogen collisions, the equilibrium equation for the three-level system (^P2, ^Pi, and ^Pq states of Ol) 
can be written as 



where /cqi and ko2 are the rate coefficients of the 0-H collisions, and Aio = 8.95 x 10 ^ s ^ is the 
spontaneous transition probability (all other coUisional and photon rates are negligibly small in this case). 

The rate coefficients for the downward transitions kio and feo, induced by 0-H collisions in the 
temperature range from 50 K to 1000 K, are given by Launay & Roueff (1977). The corresponding excitation 
rates for upward transitions can be calculated from the principle of detailed balance which implies 



where 52 = 1, 5i = 3 and go — 5 are the degeneracies of states "^Pq, '^Pi and ■^P2, respectively. AEqi = 228 K 
and AE02 = 326 K are the energy differences between the ground state and the corresponding fine-structure 
levels. 

In the observed spectrum of QSO 0000-2620 the most stringent limit on the population of the ^Pi 
state may be set by the Oi*A990.2043 A transition (see Figure 9). Using the 6-parameter value of 
10 km s~^, we estimate N{Oi*) < 1.9 x 10^^ cm~^, at 1 u confidence level. Then, with the observed 
ratio iV(0i*)/7V(0i) < 7.2 x 10"^, and assuming T ~ 10^ K (see Paper II), we find from Launay & 
Roueff kiQ = 9.22 X 10~" cm^ s'^ and fc2o = 3.80 x IQ-" cm'' or fcoi = 4.4 x 10"" cm-' s"^ and 
ko2 = 0.55 X 10"^^ cm-^ s"^, and get nni < 130 cm"''. (Note that for the chosen temperature and for 
nc/nm — 10"^ the rate coefficient for collisions with electrons equals to 3 x 10~^^ cm-^ if the collision 
strength is fioi = 3.35 x 10"^ according to Bell et al. 1998). 

Since O i* traces neutral phase regions along the line of sight, the above result may be considered as 
the upper limit for the whole Hi cloud. If this is the case, then the physical size of the z^hs = 3.3901 cloud 
is i > 7 pc. If the lower limit is nm > 1 cm~^, which is typical for the Galactic diffuse Hi clouds, then 
L < 850 pc. Comparable sizes have been also placed by Giardino & Favata (2000) for the second damped 
system at 2;abs= 3.054 towards QSO 0000-2620 from the detection of the Cii and Cii* absorptions. 



NjOi*) _ na (fcpi + fco2) 
iV(Oi) ^ Aio 



(1) 




(2) 
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4. Conclusions 

UVES observations of the quasar Q 0000-2620 gave as a result the first detection of the AriA1066.660 
A and PiiA963.801 A transitions originated in the damped Lya system at Za,hs= 3.3901. The derived 
abundances are [Ar/H] = —1.91 ± 0.09 and [P/H] = —2.31 ± 0.10 with most errors originated in the neutral 
hydrogen column density. 

The phosphorus abundance provided here is one of the few available for this element in any astronomical 
site and unique in material poor of metals. The low abundance of the non-refractory element P confirms 
the low mctallicity level of the galaxy which is therefore in the early stages of its chemical evolution. The 
P abundance is slightly below the iron-peak elements and when compared to its close nuclei S and Si, P 
shows mild evidence of an odd-even effect, hut less extreme than what predicted by theoretical models by 
Woosley & Weaver (1995) or Limongi et al. (2000). 

Ar abundance is remarkably similar to that of oxygen, sulphur and silicon, [0/H] = — 1.86 ± 0.15, [S/H] 
= — 1.98± 0.09, and [Si/H] = —1.91 ± 0.08. Ar has a large photoionization cross-section and the similarity 
of its abundance with the other a-chain elements implies the absence of significant photoionization either 
by external or internal sources. 

Log (Ar/0) and log(Ar/S) ratios are almost the same as those found in the blue compact galaxies and 
in present day galaxies where O is more abundant by one or two orders of magnitude. This strengthens the 
constancy of Ar/0 and Ar/S ratios and lends support to the existence of a universal IMF. 

In the DLA under study the relative abundances of a-chain and iron-peak elements, as measured by an 
unprecedented number of indicators, show only a mild overabundance at variance with the pattern of the 
old stellar population of the Galaxy. Since this DLA is a dust-free system, its chemical composition and the 

relative abundances of the observed species may be used to constrain model parameters in the theoretical 
calculations of stellar nucleosynthesis in the metal-poor gas. 

Measurements of the volatile and nonvolatile element abundances have shown that none of the known 
processes of dust formation or destruction is effective in the Zahs = 3.3901 galaxy. 

Finally, the upper limit to the population of the first fine-structure level in the ground state of O I is 
used to infer a lower limit to the physical dimension of the DLA system, which is placed at L > 7 pc. 
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Fig. 1.— NiiiA1709.600 A (top panel), Feii A1112.048 A (middle panel) and P ii A963.801 A (bottom panel) 
lines associated with the z^hs = 3.3901 DLA system towards QSO 0000-2620 (dots and la error bars). The 
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lines are aligned taking zero radial velocity in correspondence to Zahs = 3.390127. Smooth lines are the 
synthetic spectra obtained from the fit (see text). 
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Fig. 2.— A portion of the normalized spectrum of QSO 0000-2620 in the region of the Ar I A1048.220 A 
and AriA1066.660 A lines (Aobs = 4601.8 A and 4682.8 A, respectively) associated with the Zahs= 3.3901 
DLA. All the highest windows in the forest have been used for continuum tracing as well as other windows 
outside the portion displayed in a spectral range spanning about 1000 A. 
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Fig. 3.— ( a) An A1066.660 A line associated with the z^bs = 3.3901 DLA system towards QSO 0000-2620 
(dots and error bars). The line is aligned taking zero radial velocity in correspondence to Zabs = 3.390127. 

Smooth Uncs arc the synthetic spectra obtained from the fit as described in the text, (b) A consistent fit 
of the red wing of the AriA1048.220 A fine (the dotted curve) based on the best fit of the AriA1066 A line 
shown in panel a by the dotted curve as welL 
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Fig. 4. — Confidence regions in the Xmin ~ ^ plane calculated from the simultaneous fit of the metal absorption 
lines located redwards the Lya emission : FeiiA1611, NiiiA1709, SiiiA1808, ZniiA2026, CriiA2056, and 
CriiA2062. The parabola vertex corresponds to the best value of 6 = 10 km s~^ . 
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Fig. 5. — Results of the fitting procedure for the PiiA963 A hne. (a) The Xmin ^ ^ function of the 
column density of the singly ionized phosphorus calculated from the simultaneous fit of PllA963 and 
metals FeiiA1611, NiiiA1709, SiiiA1808, ZniiA2026, CriiA2056, CriiA2062. The 68.3% and 95.4% 
confidence levels are marked (the lower and upper horizontal dotted lines, respectively). The parabola 
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vertex corresponds to the best value of 7V(Pii) = 4.29 x 10^^ cm~^. ( b) The AC/C variations of the local 
continuum level. The grey area restrics the la deviations of AC/C in accordance to panel a), (c) The 
corresponding values of the Doppler ^-parameter. 
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Fig. 6. — As Fig. 5 but for the Ar i A1066 A line. The parabola vertex in panel a corresponds to the best 
value of TV (An) = 1.044 x lO" cm-^. 
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Fig. 7. — As Fig. 5 but for Fe ii A1112 A line and the following set of the metal absorption lines from the red 
portion of the quasar spectrum : NiiiA1709, SiiiA1808, ZniiA2026, CriiA2056, CriiA2062. The parabola 
vertex in panel a corresponds to the best value of A''(FeiiA1112) = 7.095 x 10^^ cm~^. 
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Fig. 8. — The theoretical nucleosynthetic yields (the shaded area shows ±0.2 dex uncertainty range, see text 
for more details) for the metallicity [Fe/H] = - 2.0 in the standard model of the galactic chemical evolution 
by Timmes et al. (1995). The dots with 3a error bars mark the measured relative abundances [X/Fe] at 
z^t,^ = 3.3901 towards QSO 0000-2620. 
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Fig. 9.— Oi*A990.2043 A (transition ^Pi ^T>2) expected position in the 0abs= 3.3901 DLA system 
towards QSO 0000-2620. The radial velocity scale is the same as in Figs. 1 and 3. The dotted curve is the 
convolved with the spectrograph function synthetic spectrum constraining the total O I* column density for 
the fixed value of 6 = 10 km s~^ . 



Table 1. Column densities and abundances in the dust-free pattern at Zahs = 3.3901 
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-1.91 ±0.09* 


iCr 


13.10 ± 0.03s 


-8.31 ±0.085 


-6.32 ±0.03 


-1.99 ± 0.09 <^ 


iFe 


14.87 ±0.03 


-6.54 ± 0.085 


-4.50 ±0.01 


-2.04 ±0.09*^ 


ilNi 


13.39 ± 0.03 


-8.02 ±0.085 


-5.75 ±0.01 


-2.27±0.09'^ 




12.01 ±0.05 


-9.40 ± 0.094 


-7.33 ±0.04 


-2.07±0.10« 


NOTES - 


(a) N in cm~^ 


, (b) data from Grevesse et al. 


(1996) except 



for Ar where the weighted average values from Sofia & Jenkins (1998) are 
used, (c) errors in [X/H] include errors in column densities and in solar 
abundances, (d) from Lu et al. (1996), (e) from Molaro et al. (2000), (f) 
this work, (g) mean values. 



Table 2. The [0,Si,S,Ar/Cr,Fe,Ni,Zn] ratios in the ^abs = 3.3901 system 





Cr 


Fe 


Ni 


Zn 





0.13 ±0.10 


0.18±0.10 


0.41 ±0.10 


0.21 ±0.11 


Si 


0.08 ±0.03 


0.13 ±0.04 


0.36 ±0.04 


0.16 ±0.05 


S 


0.01 ±0.04 


0.06 ±0.04 


0.29 ±0.04 


0.09 ±0.06 


Ar 


0.08 ±0.04 


0.13 ±0.06 


0.36 ±0.05 


0.16 ±0.06 



